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We show that photoassociation of fermionic atoms into bosonic molecules inside an optical lattice 
can be described using a Jaynes-Cummings Hamiltonian with a nonlinear detuning. Using this 
equivalence to the Jaynes-Cummings dynamics, we show how one can construct a micromaser for 
the molecular field in each lattice site. 

PACS numbers: 03.75.Lm, 42.50.Pq, 03.75.Ss, 42.50.Ar 
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The recent experiment by Greiner et al. that produced 
a Mott-insulator transition in ultracold ^^Rb confined in 
an optical lattice has stimulated a great deal of inter- 
est in strongly correlated systems of ultracold atoms in 
optical lattices same time, there has been a 

growing interest in the possibility of forming a molecular 
Bose-Einstein condensate (BEC) via either photoassoci- 
ation 0, 111 or a Feshbach resonance. To date, experi- 
ments have demonstrated the formation of large numbers 
of molecules via photoassociation of an *^Rb BEC and 
using a Feshbach resonance applied to an ®^Rb BEC 
and a degenerate Fermi gas of |||. 

These two lines of research have recently converged in 
theoretical proposals to photoassociate bosonic atoms in 
the Mott-insulator state into molecules 0,0,0,0,^3- 
Photoassociation in optical lattices offers distinct advan- 
tages over quasi-homogenous systems. In particular, one 
can obtain occupation numbers per lattice site that are of 
order unity, thereby minimizing the losses due to inelas- 
tic collisions. Secondly, the center-of-mass states of the 
atoms and molecules are restricted to the lowest energy 
Wannier state at each lattice site, thereby avoiding the 
difficulties associated with a multimode problem [lj| . 

In this letter we extend this work to the case of 
fermionic atoms. Our central result is that as a con- 
sequence of Pauli statistics, this problem reduces math- 
ematically within each lattice site to a situation almost 
identical to that of the familiar micromaser HKIi of 
quantum optics, the only significant difference resulting 
from the presence of an intrinsic nonlinearity in the case 
of matter-wave fields. It is well established that the ra- 
diation field generated in a micromaser exhibits a num- 
ber of quantum mechanical features that are absent in 
normal lasers and masers. For example, niicromasers ex- 
hibit strong sub-Poissonian statistics [13 and even the 
ability to dynamically generate Fock states 0|. Unlike 
lasers, which approach a classical coherent state far above 
threshold, micromasers undergo multiple phase transi- 
tions characterized by sharp changes in the photon statis- 
tics from sub-Poissonian to super- Poissonian [l^ . These 
unique features find their origin in the "granular" nature 
of the radiation field, combined with the lack of noise, 
such as spontaneous emission, which washes out the ef- 
fects of the quantum-mechanical atom-field dynamics in 



conventional lasers. 

The close analogy between the situation at hand and 
the micromaser therefore opens up exciting novel avenues 
to generate non-classical states of molecular fields and to 
study new types of dynamical quantum phase transitions 
in matter wave systems. More generally, it permits the 
extension of the numerous, well-established applications 
of cavity QED, to matter-wave optics. Of particular in- 
terest in this context is the fact that optical lattices au- 
tomatically lead to the realization of an array of individ- 
ually addressable molecular micromasers. 

We proceed by giving a broad overview of the theory 
of molecular micromasers. We discuss under which con- 
ditions the problem of fermionic photoassociation in a 
lattice can be mapped to the micromaser problem, and 
then present some of the important physical implications 
of this result. 

We consider a 2-D optical lattice in the a;?/-plane that 
confines fermionic atoms of mass m/ in two hyperfine 
states |1) and |2). The atoms are coupled to a molecular 
boson of mass mf, = 2mf via a two-photon stimulated 
Raman transition jj, |y] . The optical lattice potential 
Vf seen by the atoms is V/(r) — ^^^^ y V^^^ cos^ k^S^ + 



^,vz z'^ where V^'^'' > and we have included a trans- 
verse confining potential along the z-axis. A similar po- 
tential, with / ^ TO, confines the molecules. We assume 
that the system is at zero temperature and that the fill- 
ing factor — the number of ferniions of each type that 
occupy every lattice site — is < 1 at all times. Hence 
the fermions only occupy the lowest Bloch band of the 
lattice. For deep lattices, for which the energy separa- 
tion between the first and second Bloch bands is much 
greater than the atom-molecule interaction energy, the 
molecular state formed by photoassociation has likewise 
a center-of-mass wave function in the first Bloch band 

Si- 

The Hamiltonian for the atom-molecule system in the 
first Bloch band of the lattice is H = J2i {Hqi + Hn 
Ht where. 



Hqi = h{LUb + 6)hbi + hujf{nii + n2i) 

+h^Ubnbi{nbi - 1) + TiUxfibiifiii + n2i) -f TiUfniih2i 



2 



Hh = rix{t)b\ci.,C2., + ff.c, (1) 

and Ht accounts for tunnelling between lattice sites. 

Here, hi is the bosonic annihilation operator for 
molecules in the Wannier state of the lowest Bloch band 
centered at the lattice site i, Wb(r — r^), and similarly, Co-.i 
is the fermionic annihilation operator for state \a) with 
the Wannier wave function w/(r — r^). The correspond- 
ing number operators fiu — b\bi and ri^i — cl^Ca-^i have 
eigenvalues rihi and no-i, respectively. The terms propor- 
tional to Ub, Ux, and Uf describe the on-site two-body in- 
teractions between molecules, atoms and molecules, and 
atoms, respectively. 

The interaction Hamiltonian, Hn, describes the con- 
version of atoms into ground-state molecules via two- 
photon stimulated Raman photoassociation. x{i) is pro- 
portional to the far off-rcsonant two-photon Rabi fre- 
quency associated with two nearly co-propagating lasers 
with frequencies u)i and W2 0- Note also that the 
matter fields have been written in an interaction rep- 
resentation in which the molecular field is rotating at 
the frequency difference of the two lasers, so that 5 ~ 
Vm — {vi + V2) — (oJi — UJ2) where hv^- and hv^ are the in- 
ternal energies of the atoms and molecules, respectively. 

We assume in the following that the lattice is suffi- 
ciently deep that atomic tunnelling can be ignored dur- 
ing the time intervals r when the photoassociation lasers 
are on. We further assume that the molecules are in the 
Mott-Insulator state, for which tunnelling is completely 
suppressed for the duration of the experiment [3, l2C| . 
This second assumption is essential in the realization of 
a molecular micromaser, but we emphasize that in the 
open system that we consider in the following, it does 
not imply that the molecules at a given site are con- 
strained to be in a Fock state. Consistently with these 
assumptions, we ignore Ht altogether in the following. 

In the absence of intersite coupling, we need only con- 
sider the Hamiltonian h = Hoi + Hn at a single lattice 
site. In what follows, we drop the lattice site label for 
notational clarity. We then proceed by introducing the 
mapping ^21], 

C- — C1C2, (T^ — — clc\ 

<^z = c\ci + c\c2 - 1 (2) 

where a+ — \e){g\ and (t_ = |.9)(e| are the raising and 
lowering operators for a fictitious two-state system, and 
Uz = |e)(e| — \g){g\ is the population difference between 
its upper and lower states. Here, |e) = C2c||0) and 
\g) — |0). We note that this mapping only holds if ci and 
C2 are fermionic operators, hence, our subsequent discus- 
sion would not hold for bosonic atoms. The mapping (0) 
allows us to reexpress h exactly as, 

h = h{uJb + Ux)nb + h{ujf + Uxnb)crz 

+ n{xitpa-+x*it)ba+) + ^Ubhb{nb-l) (3) 



where we have dropped constant terms and made the 
redefinitions cob + S ^ ujb and ujf + Uf /2 ^ uif. 

In the limit Ub,Ux — > 0, this Hamiltonian reduces to 
the Jaynes-Cummings model of interaction between a 
quantized, single-mode electromagnetic field and a two- 
level atom. This model is a cornerstone of quantum op- 
tics . Because it is exactly solvable, it permits the un- 
derstanding of detailed aspects of the dynamics of light- 
matter interaction. The development of ultra-high Q cav- 
ities and cavity QED has led to the experimental realiza- 
tion of systems that closely approach the ideal situation 
of the Jaynes-Cummings model. Of particular relevance 
in the present context is the micromaser, which consists 
of an ultra high-Q microwave cavity in which two-level 
atoms are injected one at a time. 

The dynamics of the micromaser are governed by three 
mechanisms: the injection of a sequence of individual 
atoms from a very dilute atomic beam inside the mi- 
crowave cavity, the Jaynes-Cummings interaction be- 
tween these atoms and the cavity mode, and cavity dis- 
sipation, which can normally be neglected while an atom 
transits through the cavity, but that dominates the field 
dynamics in the intervals between atoms. A similar sit- 
uation can be achieved in the present system: all that 
is required is to inject a sequence of pairs of fermionic 
atoms inside the lattice well, turn on the photoassocia- 
tion lasers for some time interval r to introduce a Jaynes- 
Cummings-like dynamics — where the electromagnetic 
field mode is replaced by the molecular field — and fi- 
nally turn off these fields for a time T and let dissipation 
take over. This sequence is then repeated to build up 
the molecular field. The close analogy with the micro- 
maser makes it clear that this "molecular micromaser" 
can generate strongly nonclassical molecular fields, and, 
more generally, opens the way to "cavity atom optics." 

We now discuss the three key mechanisms that gov- 
ern the "molecular micromaser" dynamics, starting 
with the Jaynes-Cummings-like evolution. For x = 
const, equation can be solved just like the Jaynes- 
Cummings Hamiltonian within the two-state manifold 
{|e, Ub), \g, Ub -f 1)}. The dynamics is then in the form of 
quantized Rabi oscillations between the two states at the 
frequency 

= \j [2ujf -LUb + {2Ux - Ub) Ubf + 4|x|2(nb -t- 1). 

(4) 

In particular, if the system is initially prepared in the 
state \e,nb), the probabilities for the system to be in the 
two states of the manifold after a time r are |ce,n^ (t) | = 
1 - C„,+i(t) and |cg,„,+i(T)|^ = C„,+i(t) where 

C„(r).^sin^ Q7^„_,r). 

Note that unless Ub — 2Ux, the detuning in TZn depends 
on the number of molecules present and thus resonant 
Rabi oscillations are only possible for a single manifold. 
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The pump mechanism consists of continuously inject- 
ing pairs of fermionic atoms into each lattice site dur- 
ing the interval T when x = 0. This can be realized 
in a 2-D lattice by a beam of fermions incident on the 
lattice along the z-axis in internal states 1/^) that have 
polarizabilities opposite to those of \a = 1,2). This re- 
sults in them seeing a repulsive instead of an attractive 
potential centered at z = 0. The energy of the inci- 
dent beam is taken to be much less than the barrier 
height, so that tunnelling is negligible. Two lasers di- 
rected along the lattice can then be used to stimulate 
Raman transitions from the untrapped to the trapped 
states, l/o-) — !■ \a). The coupling constants are then given 
by iKk = / d^re*'''^ri2(r)ic}(r)/2\/T^ where Q2 is the two- 
photon Rabi frequency and V the normalization volume 
for the plane wave states of the incident beam 0. 

For a broadband continuum of reservoir states, one can 
derive a Born-Markov master equation for the density 
matrix of the trapped fermions, 



to both loss mechanisms is 



dpit) 



dt 



^ -9 (" t^- 



2ct 



+ (1 - »^) [4c<tP - 2c„p4 + pclca] ) (5) 



where n is the number of incident fermions with en- 
ergy fik^/2mf = LUf and the "pumping rate" T — 
2tt\k{u! f)]"^ D{uj f) with D{uj) being the density of states 
for the plane waves |22l |. In order to pump the lat- 
tice site with fermions we take h — I. Since the two 
fermionic species evolve independently in Eq. (jSJ, wc 
can express their density matrix as the tensor product of 
the density operator for each species, p'^^\t)®p'''^\t). Eq. 
(O then has the solution p\^l{t) = cxp(— r/;)pQ'^Q (0) and 
pS(^) = 1 - p["l{Q)eM'Tt) where p'^^ = (0|p('^)|0) 
and p^^l = (0|c^p(^)4|0). For Tt > 1, p^^l{t) = 1, i. 

e. the fermions are in the state |e) — C2c||0) with unit 
probability. 

Finally, the damping of the molecular mode is de- 
scribed via a phenomenological master equation that 
includes the dominant loss mechanisms. They include 
in particular three-body inelastic collisions between a 
molecule and two atoms in which the atoms form a dimer 
with the resultant loss of the atoms and molecule. This 
occurs at a rate 73 ^ ha^^ / m f{l\i'^ (.[J^ l^^^Y where a/ is 
the atomic scattering length and l'"^^ the characteristic 
spatial extent of w/ (r) along the three axes jU |3| ■ An- 
other loss mechanism is spontaneous emission from the 
molecular excited state from which the optical lattice and 
photoassociation beams are detuned. When x(0 = 0; 
these losses are due solely to scattering of lattice pho- 
tons. This gives a rate 71 ~ 7eW^''V4|A| for a blue de- 
tuned lattice, using a harmonic approximation for the 
optical potential at the lattice nodes. Here 7e is the ex- 
cited state linewidth and A the detuning from the ex- 
cited state. The contribution to the master equation due 



dp{t) 



dt 



loss 



--^[h^hp^ 



6p6t] [B^Bp - BpBl] +H.C. 

(6) 

where B — bciC2 = ba^ 24]. For molecules created in 
their rotational-vibrational ground state, we can ignore 
losses due to rotationally or vibrationally inelastic colli- 
sions between molecules, or between atoms and molecules 
la]. _ 

In the limit of strong pumping, F 3> 71,3, the atoms 
reach their steady state before the state of the molecules 
has noticeably changed. We can then replace Eq. © 
with a coarse-grained master equation valid for time in- 
tervals much larger than F"^. To do this we substitute 
the steady-state values of Eq. © into the density oper- 
ator, p{t) = p^^\t) (g) |e)(e| where p^^^ is the molecular 
density operator, and then trace over the states of the 
atoms in Eq. ©. This yields finally 



1 



dt 



b^bp'-''^ - 2bp^^'^b'^ + pC-^fetfo = 



where 7 = 71 -I- 73. We note that for aj '-^ lOOao and 
TO/ ~ 10a. TO. M. one has 73"^ ^ 10s while 7]"^ ~ Is for a 
detuning of 10000 hue widths and w^') - 27r x lO^s'^ 

Just like in micromaser theory, the weak damping of 
the molecular mode allows us to assume that during the 
intervals r of photoassociation the evolution is unitary 
and given by h. Since r ~ l/VrT+Tlxj for complete 
Rabi oscillations to occur in the presence of n molecules, 
this implies that the photoassociation fields must satisfy 
Vn+l|x| 3> 71,3. The resulting molecular density oper- 
ator is then 



pW(t/+T) =Tr, 



— ih 



llT _ 



F{r)[pm (7) 



where ti = 1{T + r) for I = 0,1, 2, .... In order to deter- 
mine F{t) explicitly, we recall that at the end of each 
interval T every lattice site is occupied by a pair of 
fermions with unit probability for FT 1. Since the 
fermions are uncorrelated with the state of the molecules 
already present, the reduced density matrix at + t is 
then p('')(i;-fr) = I]„_Pn(ti +r)|n) (n| where p„ (t; -f r) = 
[l-Cn+iir)]Pniti)+Cn{T)pn-iiti). Here p„ = (n|pW|n) 
and we have assumed that p^''^ is initially diagonal in the 
number state basis. 

The formal solution for the molecular density operator 
is then = exp(£r)i^(r)p('')(t/). We note that T 

is controlled by the timing of the photoassociation beams, 
in contrast to conventional micromasers where T corre- 
sponds to the arrival times of the atoms and is therefore 
normally a random variable. 

The steady-state molecular field is given by the return 
map condition p^**) = = p^''\ti). If the interval 

T is such that F~^ ^ T ^ 7^^, we can expand the 
exponential in the formal solution to lowest order in C. 
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FIG. 1: Average number of molecules, {ni,}, versus Q for 
iVe^ = 25, 77 = 0, and = O(solid line), 0.05 (dashed dot), 

0. 15 (dotted line). 

The steady-state condition is then C[p^''^ = r~^(l — 
F{t))p^''\ which has the same form as the steady state 
solution of the micromaser if one identifies with the 
rate at which atoms are injected into the cavity |l5|. The 
steady state solution in the number-state basis is, 

P--Pof[^ (8) 

1=1 ' 

where pn is the probability of having n molecules and po is 
determined by the conservation of probability, J^nPn = 

1. Eq. |(Hl has the same form as the zero-temperature 
micromaser photon statistics, except that C;(r) contains 
an Z-dependent detuning in TZi resulting from the two- 
body interactions. 

Figures 1 through 3 highlight important features of the 
molecule statistics as a function of Nex = 1/ jT, the num- 
ber of pumping cycles per lifetime of the molecules, and 
of the dimensionless strength for the nonlinear detuning 
f3 = {2Ux — Ub)/2\x\- The micromaser pump parameter 
is given by Q — {N^xY^'^X'''- Figure 1 shows the average 
number of molecules at each site, as a function of 
for the linear resonance condition, 77 = (2ajy — ti;f,)/2|x| — 
0. We see that the lasing threshold is not affected by 
the nonlinear detuning and occurs at « 1. However, 
above threshold {hb) is strongly suppressed as (3 is in- 
creased. This is because the effective detuning between 
the fermion atom pairs and molecules increases with in- 
creasing molecule number due to the self-phase modula- 
tion of the molecules, Ubnh{nb ~ 1), and the AC Stark 
shift of the atoms, UxUhdz- For /3 > 0(< 0) this can be 
compensated for by choosing a finite negative (positive) 
detuning, 77 < 0(> 0), as shown in Fig. 2. 

Figure 3 shows the Mandel Q parameter, Q — 
(("-fc) ~ (^&)^) /("-f)) ~ 1- Just above threshold, the 
molecule distribution is strongly super-Poissonian (Q > 
0) and then become sub-Poissonian {Q < 0) until 8 « 27r 
after which Q shows very sharp oscillations. We note 
that the number fluctuations decrease with increasing 
\P\ and show smaller super-Poissonian peaks. The sharp 
resonances in Q and (ni,) for large Q are attributable 
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FIG. 2: {m,) versus 6 for iV,^ = 25 and f3 = ^005 for 77 = 
(solid line), -1-1 (dashed dot), —1 (dotted line). 
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FIG. 3: Q versus 9 for iV,^ = 25, 77 = 0, and (3^ = (solid 
line), 0.05 (dashed dot), 0.15 (dotted line). 



to downward trapping states [ijj- Unlike the conven- 
tional micromaser, there is no thermal noise present in 
the damping of the molecular field. If we replace /^[p*-^'] 
with the master equation for coupling to a finite tem- 
perature bath with a thermal occupation of hh, then the 
rapid fluctuations in Q are averaged out. 

In conclusion, we have examined the photoassociation 
of fermionic atoms into bosonic molecules and shown that 
the dynamics can be described by an effective Jaynes- 
Cummings Hamiltonian with a nonlinear detuning re- 
sulting from the two body elastic interactions of the 
molecules. We have also presented the theory of a molec- 
ular micromaser. 
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